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Abstract 

The standard electroweak model is extended by means of a second Brout- 
Englert-Higgs-doublet. The symmetry breaking potential is chosen is such a 
way that (i) the Lagrangian possesses a custodial symmetry, (ii) a static, spher- 
ically symmetric ansatz of the bosonic fields consistently reduces the Euler- 
Lagrange equations to a set of differential equations. The potential involves, 
in particular, products of fields of the two doublets, with a coupling constant 
A3. Static, finite energy solutions of the classical equations are constructed. 
The regular, non-trivial solutions having the lowest classical energy can be of 
two types : sphaleron or bisphaleron, according to the coupling constants. A 
special emphasis is put to the bifurcation between these two types of solutions 
which is analyzed in function of the different constants of the model, namely 
of A3. 
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1 Introduction 



It is known for a long time [1] that baryon and lepton numbers are not strictly con- 
served in the standard model of electroweak interactions (see [2], [3] for reviews). 
Baryon number violating processes [4] involve the crossing through an energy barrier 
separating topologically inequivalent vacua of the underlying gauge theory. Remar- 
quably, this energy barrier is high but finite. It corresponds to a static, regular 
solution of the classical equations of motion : the sphaleron [5] . The sphaleron was 
first constructed in the case 9w = {6w denotes the Weinberg angle) where a con- 
sistent spherically symmetric ansatz [6] transforms the Euler-Lagrange equations of 
the theory into differential equations. The Klinkhamer-Manton (KM) sphaleron is, 
however, not the minimal energy barrier when the mass of the Brout-Englert-Higgs- 
boson (BEH-boson) exceeds some critical value. Indeed, for Mh » Mw [7, 8, 9] 
another branch of solutions exists which bifurcate from the sphaleron branch for 
Mh ~ 12Mw- The new solutions have a lower energy than the sphaleron and since 
they appear by pairs connected each other by the parity operator, they were called 
bishpaleron. Nowadays the possibility that bisphalerons constitute the energy bar- 
rier allowing for baryon number violating process is ruled out in the minimal (one 
doublet) electroweak model by the perturbative upper limit of the BEH-field and the 
sphaleron-bisphaleron bifurcation remains a curiosity of the classical equations. 

However several extensions of the minimal Wcinbcrg-Salam model arc curently 
under investigation as alternative candidates for the description of electroweak inter- 
actions (see e.g. [10, 11]). Among these various extensions, the ones incorporating 
more than one multiplet of BEH-boson play a central role. For instance, the minimal 
supersymmetric electroweak model, considered for many theoretical reasons, involves 
two BEH-doublets. These extended models lead generally to involved classical equa- 
tions where the generalisations of sphaleron and bisphaleron can be looked for, as well 
as eventual other type of solutions of soliton type. It particular, it is challenging to 
study the domain of parameters for which bisphaleron exist and to see if it intersect 
the domain of physically acceptable parameters. This question was adressed early in 
[12] and in [13]. The potential used these papers does not involve a direct coupling 
between the doublets. Here we will extend the potential choosen in these two pa- 
pers by mean of a supplementary interaction between the two BEH-doublets. The 
infiuence of the new term on the sphaleron-bisphaleron bifuraction will be studied in 
details. 

To be complete let us mention that the classical equations of the two-doublets- 
extended standard model were also investigated in [14, 15] with even more general 
potential but, to our knowledge, these authors did not put the emphasis on bifurcation 
between the two types of lowest energy solutions. 

In Sect. 2, wc present the model, the notations and the physical parameters. The 
spherically symmetric ansatz, the equations and boundary conditions are given in 
Sect. 3; the numerical solutions are then discussed in Sect. 4. 
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2 The model 



The lagrangian that we consider in this paper reads 

£ = -If^^F^^''^ + + (L>^$(2))t(L>''$(2)) - $(2)) (1) 

Where $(2) denote the two BEH-doublets and the standard definitions are used 
for the covariant derivative and gauge field strengths : 

C = d,V: - d,V^ + ge'^'XK (2) 

D,^ = {d, - '-gr'^V^)^ (3) 

(the hmit 9w — 0, i.e. g' = 0, is used throughout the paper). 

The most general gauge invariant potential constructed with two BEH-doublets is 
presented namely in [11], it depends on nine constants. Here we consider the family 
of potentials of the form 

2 2 22 

V^($(l),$(2)) = Ai(<l>{,)$(l)-^)2 + A2($j2)$(2)-|)' + A3($fi)$(l)-|)(^^^^ 

(4) 

depending on five parameters. The terms directly coupling the two doublets is 

parametrized by the constant A3 and one of the main characteristic about this poten- 
tial resides in the fact that it imposes a symmetry breaking mecanism to each of the 
BEH-doublcts. The case A3 = is studied at length in [12, 13]. 

The lagrangian (1) is off course invariant under SU(2) gauge transformations but 
it further possesses a larger global symmetry under SU(2) x SU(2) x SU(2). In fact, 
the part of the lagrangian (1) involving the scalar fields can be written in terms of 
2x2 matrices defined by 

When written in terms of the matrices Mi and M2, the lagrangian (1) becomes 
manifestly invariant under the transformation 

V^ = AV^A^ , M[ = AMiB , M^^ = AM2C (6) 

with A,B,C e SU(2); this the custodial symmetry. The double symmetry breaking 

mechanism imposed by the potential (4) leads to a mass for two of the three 
gauge vector bosons and, namely, to two neutral BEH-particles with masses M^, Mh- 
In terms of the parameters of the Lagrangian, these masses are given by [11, 16] 

Mw = yv! + v'2 , Ml, = l[Ai + ^2 ± - ^2)2 + 452] (7) 
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with 

A, = 2vfi\i + \s) , A2 = 2vl{X2 + X3) , B = 2XsViV2 (8) 
For later convenience we also define 



tan/3 = — , phj^ = —— , ep = 4— , p= 1,2,3 (9) 
vi Mw 5 



Note that the mass ratio pi^2 used in [13] are related to pH,h by Ph — niax{pi,p2}, 

Ph = min{pi,p2}- For physical reasons, we consider only > 0, f2 > so that 
< j3 < 71 / 2. Interestingly, the parameter 63 can be negative but cannot take 
arbitrary values. The following relations are usefuU to determine the physical region : 

ei cos^ [3 + 62 sin^ /3 = ^pjj + pi - 63 

eiCOs'/3-e2sin2/3 = (pjj - plf - 4ei sin\2(5) - £3 cos(2/?) (10) 
The physical domain is then determined by the conditions 

'"'"''^ <^3<^^. > ^<pI<p'h (11) 



2 sin(2/3) " ^ " 2 sin(2/?) 
Physical domains are presented on Fig.l for the case = 3 for /3 = 7r/4 and P — Sn/S. 

3 Spherical symmetry 

In order to construct classical solutions of the Lagrangian (1), we perform a spherically 
symmetric ansatz for the fields. With the notations of [17], it reads 

Vo" = 

..a . . fsir) _ fcir) ^ ^ 
Vi = eaijrj H {Sia - riTa) H r^r^ 

gf gf gf 



0(1) = ^[i/(r) + i(f.a)X(r)](^° j 



0(2) = ^[^(r) + i(f.a)X(r)](^ J j (12) 

Where /a, /b, -Fc, H, K, H, K are real radial functions. It can be shown that the 
above ansatz transforms the Eulcr- Lagrange equations into a set of coupled differential 
equations. The custodial symmetry has been used to set the two doublets parallel 
to each other asymptotically. The condition Vq = results from a gauge fixing. In 
fact, the spherically symmetric ansatz leaves a residual gauge symmetry which can 
be exploited to eliminate one of the seven radial functions [7, 8, 17]. Here we will 
adopt the radial gauge XjV-^ — which implies fc — 0. 
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The effective one-dimensional energy density can be obtained after some algebra 



E ^ ^ redx 

aw Jo 

^ (13) 

aw 



2 

with aw = 4^ and 



+ cos' P[{H{U - 1) + KfsY + {K{fA + 1) - Hfsf + 2x\{H'f + {K'f) 

+ {H{U - 1) + KfB? + {K{U + 1) - Hfsf + 2x\{H'f + {K'f)] 

+ cos^ /3[eix^(//^ + - 1)2 + 62X^(^2 + - i//3)2 

+ esx''{H' + K''-l)(H' + K''-tg''P)] (14) 

The dimensionless variable x = Mwt is used and the prime means derivative with 
respect to x. The equations to solve can then be obtained by varying the functional 
(14) with respect to the six radial functions. Remark that in the case V2 = 0, A2 = 
A3 = the equations for H, K decouple and these functions can be set consistently 
to zero; the remaning system just corresponds to the one of the IDSM. 

The conditions of regularity of the solutions at the origin imposes namely = 
1 at X = 0, the custodial symmetry (6) can then be exploited to fix the following values 
of the radial fields at the origin 

/a(0) = 1 , /b(0) = , K{Q) = Q , k{0)=0 , H'{0) = , H'{0) = (15) 

On the other hand, the condition of finiteness of the classical energy imposes the 
following asymptotic form 

if A, fB)x=oo = (cos 2nq, sin 2nq) 

{H,K)-^=oo = (cos7r(g - /c),sin7r(g - fc)) 

{H,K)x=oo = tan/3 (cosTTQ'jSinTrgf) (16) 

for some real number q and k equal to zero or one. For later use we define g = 1/2 + 5. 



4 Discussion of the solutions 

In order to make the following discussion self contained, we first summarize the main 
features of the solutions available in the one doublet standard model (IDSM), i.e. in 
the case = A2 = A3 = with H — K — 0. 
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4.1 IDSM 



There exist at least one solution, the Klinkhamer-Manton (KM) sphaleron, for all 
values of pi [6] { pn = pi in this case). For this solution one can further set Jb = 
= by an appropriate choice of the custodial symmetry; the classical energy 
increases monotonically as a function of pi : 

Es{pi = 0) a; 3.04 , Es{pi = oo) ^ 5.41 (17) 

The KM sphaleron is always unstable but the number of its directions of instability 
increase when p increases [8, 18]. At p ^ 12.04 a couple of new solutions, the bis- 
phalerons, bifurcate from the sphaleron. The two bisphalerons (which transform into 
each other by parity) have the same energy and their energy is lower than the one of 
the KM sphaleron 

^fe,(pi = 12.04) = ^,(pi = 12.04) ^ 4.86 , EUPi = oo) ^ 5.07 (18) 

The parameter q defined in (16) is equal to 1/2 for sphaleron. For bisphaleron, the 
deviation from q = 1/2 parametrized by 6 varies from zero (at the bifurcation point) 
and S — ±0.06 (at pi = oo), respectively for the two bisphalerons. 

4.2 2DSM, case 63 = 

Solutions of the Lagrangian (1) with A3 = were first constructed in [12] and recon- 
sidered in [13] where the emphasis on the sphaleron-bisphaleron bifurcation was set. 
As in IDSM, sphaleron have fs = H = H = and seem to exist for all values of 
the parameters of the potential. The angle parameter q defined in (15) correponds to 
q = 1/2, irrespectively of the coupling constants of the potential. 

By contrast, the six radial functions corresponding to the bisphaleron are non 
trivial and fuUfil the boundary conditions (16). The parameter q depends of the 
various coupling constant although remaining close to 1/2 (e.g. 5 — 0.036 for pi — 
14,p2 = l,/9 = 0.2). 

The results of [13] show the existence of a smooth surface in the Pi, P2- /^-parameter 
space inside of which only sphaleron solution exist while sphaleron and bisphaleron 
coexist outside, the bifurcation taking place on the surface. The critical surface can 
be determined only numerically by studying a few parameters characterizing the bis- 
phaleron solutions, namely the values 6, K{(}), K{0), as functions of f3, €1^2,3- Varying 
one of these parameters and fixing the other three, a critical point is determined 
when 6, K{0), K{0) approach zero. This is done on Fig. 2 for /3 = 7r/4 ,£3 = and 
Ph — ^Ph- The critical point then corresponds to pn ~ 2.397. It is worth noticing 
that bisphaleron solutions with 5 = but K{x) 7^ 0, K(x) ^ also occur outside 
the critical surface. This is illustrated on Fig. 3 where we set pf^ + pf^ = 9 and vary 
= atan{pH / Ph)- Clearly, 6 goes to zero in the limit G — >■ 7r/4. The figures 2 and 3 
are complementary to the ones presented in [12, 13] 

The main feature of the bisphaleron solutions in the 2DSM is that the angle 
(j) — axctan{K/H) increases monotonically from (for x — 0) to tt — q (for x — 00) 
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while = aYctan{K / H) decreases from to 7r(g — 1); that is to say that these solutions 
of lowest energy have = 1 in (16). In fact, non trivial solutions with k = were 
discovered in [12] but since they have a higher energy, they are likeky less interesting 
as far as the energy barrier is concerned. 

4.3 2DSM, case 63 ^ 

The construction of solutions in the domain ph, pn, ^,€3 and the study of the critical 
hypersurface in this four-parameter space is a vast task. For definiteness, we have 
limited our investigations to the case ph — 3.0 and to two values of /3, namely (3 — tt/A 
and P = 3n/8. 

We first discuss the results for P = n/i, i.e. vl = V2. Here are a few "points" on 
the bifurcation line correponding to (3 = tt/4 : 

{ph,pH) ■■ (10,7.6) , (3.0,7.08) , (5.0,6.08) , (5.585,5.585) (19) 

Our main concern is to determine how this domain of the Ph, Ph plane evolves with 
€3. The result is illustrated by Fig. 1. where the physical domain is delimited by the 
solid lines. The sphaleron-bisphaleron bifurcation line is represented by the solid line 
with bullets and the domains where sphaleron only exist is indicated, as well as the 
domain where sphaleron and bisphaleron coexist. The critical line separating these 
two regions clearly exhibits two different behaviours on the domain of parameters 
considered : for 6.2 < pij < 8, the critical line is roughly a function increasing 
linearly with the coupling constant €3. As a consequence, the critical line intersect 
the lower line delimiting the physical domain, for instance at pn ~ 6.22 , £3 ~ —15.0. 
Clearly, for the negative values of £3, the the minimal mass of Mh (with all other 
parameters fixed) for which bisphaleron exist is lowered by the presence of a direct 
interaction between the BEH-fields. For 8.0 < pn < 10.0 the critical lines develops a 
plateau at 63 10.0 and pH,cr depends only weakly of 63. For some unknown reason, 
the numerical analysis become very difficult when the critical line appraoches the 
limit of the physical domain. 

The behavior of the critical line turns out to be completely different for /3 = 27r/8, 
in this case, the limit of the physical domain is indicated by the dashed lines and the 
critical line bu the dashed line with the triangles. In contrast to the case /3 = 7r/4 we 
see here that the critical value pH,cr decreases roughly linearly with 63. For £3 > the 
barrier turns out to be a bisphaleron for lower values of pu than in the €3 = case. 
Here, the critical line seems to cross the two lines determining the physical domain 
where it naturally terminate. 

We further studied the critical phenomenon for 37r/8 < /3 < 7r/4 and observed a 
smooth evolution of the critical lines displayed in Fig. 1. 
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5 Conclusion 



The lagrangian considered in this note leads to a tricky system of six differential 
equations with boundary conditions and depending effectively on four parameters. 
Many types of non-trivial solutions can be constructed numerically [12] but, at the 
moment, the ones with lowest energy are identified as the sphaleron or the bisphaleron, 
depending on the different coupling constants. The determination of the critical 
hypersurface of bifurcation in the space of parameters constitutes a huge task which 
can be studied only numerically. The problem is for a large part academical; however 
at the moment, the theoretical limits on the BEH-boson masses obtained in the two- 
doublets extension of the electroweak model, do not yet exclude that the energy 
barrier between topologically different vacua could be determined by a bisphaleron. 
A few years ago, it was already pointed out that the minimal mass of the neutral 
BEH-fields for the barrier to be of the bisphaleron type is considerably lower in the 
two-doublets model (without direct interaction of the doublets in the potential) than 
in the minimal, one-doublet model. The calculations reported here suggest that, if a 
custodially-invariant coupling term is supplemented to the potential, the critical mass 
PH,cr varies roughly linearly with the new coupling constant £3. The supplementary 
coupling constant cannot off course be arbitrarily large, but on the physical domain, 
the critical value pH,cr indeed dimnishes, all other masses being fixed. Off course, 
many other terms can be added in the potential [11] and it could be that bisphaleron 
solution exist for still lower masses of the neutral BEH-ficld when terms allowing for 
charged BEH-fields are considered as well. 

Acknowledgements This work started from several discussions with Michel Her- 
quet. I gratefully acknowledge him for these discussions, his interest in the topic and 
for numerous private communications about his own reseach [16]. 
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Figure 1: The physical domain corresponding to the case ph — 3.0, P — 7r/4 (resp. 
/3 = 37r/8) are represented by the sohd (resp. dashed) hne in the Ph,^3 plane. The 
line with bullets (resp. triangle) represents the sphaleron-bisphaleron bifurcation. 
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Figure 3: The evolution of the parameters 5,Ki = K{0),K2 = K{0) is shown as 
function of pn / Ph for + p| = 9. 
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